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Abstract. Significant in-band stray light (IBSL) error at so- 
lar zenith angle (SZA) values larger than 77° near sunset in 
4 SBUV/2 (Solar Backscattered Ultraviolet) instruments, on 
board the NOAA-14, 17, 18 and 19 satellites, has been char- 
acterized. The IBSL error is caused by large surface reflec- 
tion and scattering of the air-gapped depolarizer in front of 
the instrument’s monochromator aperture. The source of the 
IBSL error is direct solar illumination of instrument com- 
ponents near the aperture rather than from earth shine. The 
IBSL contamination at 273 nm can reach 40 % of earth radi- 
ance near sunset, which results in as much as a 50 % error 
in the retrieved ozone from the upper stratosphere. We have 
analyzed SBUV/2 albedo measurements on both the dayside 
and nightside to develop an empirical model for the IBSL er- 
ror. This error has been corrected in the V8. 6 SBUV/2 ozone 
retrieval. 


1 Introduction 

Stratospheric ozone shields Earth’s biosphere from solar ul- 

traviolet radiation (DeFabo, 2000), while tropospheric ozone 

as a pollutant impacts global health and economy (Selin et 

al., 2009). Atmospheric ozone plays an important role in cli- 

mate change (Fiore et al., 2002). The Solar Backscattered Ul- 
traviolet (SBUV) program for monitoring atmospheric ozone 
profiles has spanned more than three decades (Stolarski and 
Frith, 2006). The first SBUV instrument was launched on 
the Nimbus-7 satellite into a sun synchronized polar or- 
bit in 1978 (Heath et al., 1975). It has been followed by 

7 SBUV/2 instruments on NOAA satellites from 1985 to 
the present day (NOAA-9, NOAA-11, NOAA-14, NOAA- 
16, NOAA- 17, NOAA- 18 and NOAA- 19) (Frederick et al., 
1986; McPeters et al., 2013). The last 4 SBUV/2 instruments 


are still active. SBUV/2 instruments measure backscattered 
terrestrial radiance in the nadir direction at 12UV wave- 
lengths from 252 to 340 nm with a nominal spectral bandpass 
of 1 . 1 3 nm FWHM (full width at half maximum) for ozone 
monitoring. The incident solar irradiance at the top of atmo- 
sphere is also measured by SBUV/2 instruments on a regular 
basis for calibration. Measured ratios of the earth radiance 
to the solar irradiance (albedo) are used to derive ozone pro- 
files (Bhartia et al., 1996, 2013). The long-term SBUV ozone 
record provides important information to scientific research 
in climatology and meteorology, and to environmental pro- 
tection efforts (Stolarski and Frith, 2006; Terao and Logan, 
2007). 

In the SBUV/2 spectral range, the earth radiance inten- 
sity varies by more than 3 orders of magnitude. Figure 
1 shows typical spectra of radiance and albedo measured 
by an SBUV/2 instrument in its continuous spectral scan 
mode, which were taken and averaged in latitude bands from 
15° S to 15° N. The square signs indicate the 12 discrete 
SBUV/2 wavelengths for the regular ozone monitoring oper- 
ation. Light at the short wavelength channels (252-298 nm), 
mostly absorbed by ozone before reaching ground, provides 
information of ozone density in the stratosphere (Bhartia et 
al., 2013). The measurements at longer wavelengths, which 
can reach the surface, provide limited information about the 
ozone distribution at low altitude. Detecting a 2 % ozone den- 
sity change in Umkehr layers in the stratosphere requires 
the calibration of SBUV/2 albedo measurements at the short 
wavelengths to be accurate to within 1 %. With such a large 
dynamical signal range and stringent accuracy requirement, 
stray light often becomes a difficult issue in optical design 
and characterization of a spectrometer. There are 2 types of 
stray light, out of band stray light (OBSL) and in-band stray 
light (IBSL). OBSL represents contamination of the signal 
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Fig. 1. Solar backscattered ultraviolet spectrum measured by 
NOAA-16 SBUV/2. The monochromator completes a continuous 
spectral scan from 405 to 160nm in 168 s. The radiance spec- 
trum (top panel) is an average of 60 scans between the latitudes of 
15° S and 15° N, which were taken during four consecutive orbits 
once per week in November 2001. The albedo spectrum (bottom 
panel), equal to the radiance spectrum divided by the solar irradi- 
ance spectrum, is used in the ozone retrieval. Square symbols rep- 
resent SBUV/2 discrete wavelengths used for ozone measurements 
in a 32 s scan cycle. The diamond shows the location of the cloud 
cover radiometer (CCR) measurements. 


at the measured wavelength from light at other wavelengths 
due to imperfect spectral filtering. In the SBUV ozone mea- 
surements, a typical OBSL error for bright surface scenes 
can add a few percent error in radiance at the shortest wave- 
length (DeLand et al., 2012). A simplified description of the 
OBSL characterization is presented in that paper, but a de- 
tailed discussion would require substantially more space. We 
will not discuss the characterization of OBSL in this paper. 
IBSL is considered to represent imperfect spatial filtering of 
external illumination that passes through the spectral band- 
pass of a given wavelength channel. Analysis of stray light 
in a remote sensing instrument is often an important part of 
sensor characterization. When an IBSL error is present, the 
sensor’s measurements are contaminated with light that does 
not come from the research target. The signal is an additive 
error that cannot be canceled in the albedo ratio in the ozone 
retrieval. As will be described later in this paper, the IBSL 
error is also wavelength dependent as well as time depen- 
dent. Typical retrieval algorithms for atmosphere trace gasses 
could not reduce or eliminate such additive errors in their 
data products. Therefore, detection, characterization and cor- 
rection of IBSL, if significant, are necessary. We hope that 
our study of SBUV IBSL can also provide lessons for the 
design of future instruments. 

SBUV/2 instruments have an 11° (FWHM) field of view, 
which projects to approximately a 168 km x 168 km square 
on the ground in the nadir direction. The response over this 


NOAA-18 SBUV/2 Earth View Signal at 273.7 nm (2005/264) 



NOAA-16 SBUV/2 Earth View Signal at 273.6 nm (2005/264) 



Fig. 2. Comparison of SBUV/2 earth view (radiance) signals at 
273 nm between NOAA-18 (top panel) and NOAA-16 (bottom 
panel). The IBSL error in NOAA-16 data is negligible. NOAA-18 
SBUV/2 data have significant IBSL errors at SZA >77° approach- 
ing the terminator near North Pole. NOAA-16 and NOAA-18 orbits 
were in proximity at the Equator crossing local times of 14:59 and 
14:01, respectively, to minimize diurnal difference in the earth radi- 
ance. 


field of view was characterized in the SBUV/2 prelaunch cal- 
ibrations. The instrument’s sensitivity varies less than 10% 
within a 10° x 10° center area, then, decreases sharply by a 
factor of 10 -4 within 2° outside the center area. Such perfor- 
mance makes the SBUV radiance measurements well defined 
in geolocation. Prelaunch testing did not fully simulate the 
solar light direct illumination of surfaces on the spacecraft 
near the SBUV/2 instrument’s aperture, which happens at an 
incidence angle nearly perpendicular to the instrument’s opti- 
cal axis. The solar flux is 4 orders of magnitude more intense 
than the earth radiance at 252 nm, as shown by the albedo 
values in the bottom panel in Fig. 1. Figure 2 shows sample 
albedo measurements at 273 nm on the same date by NOAA- 
18 SBUV/2 (N18) in the top panel and NOAA-16 SBUV/2 
(N16) in the bottom panel. We will use the abbreviation N9, 
N1 1, N14, N16, N17, N18 and N19 for SBUV/2 instruments 
on board the NOAA-9, 11, 14, 16, 17, 18 and 19 satellites, 
respectively. Both N16 and N18 were emerging from dark- 
ness into daylight in the Southern Hemisphere and moving 
into darkness in the Northern Hemisphere in the same week 
around the fall equinox in 2005, in orbits with Equator cross- 
ing local times of 14:59 and 14:01, respectively. We selected 
two instruments with similar orbits to minimize radiance dif- 
ferences due to diurnal ozone variations. The N16 signal 
monotonically decreases in both hemispheres as solar zenith 
angle (SZA) increases and approaches zero at SZA = 95° 
and beyond. This is consistent with expectations for a short 
wavelength (A. < 300 nm) that is blocked by ozone absorp- 
tion before reaching ground. The N 1 8 Southern Hemisphere 
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DEPLOYABLE CALIBRATION 



STRAY LIGHT STOP 


Fig. 3. Diagram of SBUV/2 instrument optics. SBUV(/2) instru- 
ments consist of double Ebert monochromators for scanning ozone 
absorption wavelengths, and a bandpass-filter photometer as a 
cloud cover monitor. A depolarizer is installed in front of the 
monochromator to minimize the polarization effect with atmo- 
spheric Rayleigh scattering coupled with strong polarity with the 
grating spectrometer. Courtesy of Ball Aerospace and Technologies 
Corp. 


data are consistent with N16, but the N18 Northern Hemi- 
sphere data show a sharp step increase between SZA = 78° 
and 80° in comparison with N16. The increased signal is still 
present in nightside data until the satellite passes into Earth’s 
shadow (eclipse) at SZA =118°. This behavior suggests that 
the source of the anomalous signal is direct solar illumination 
of the SBUV/2 instrument. We found significant IBSL con- 
tamination in the earth radiance measurements at the short 
wavelengths with N14, N17, N18 and N19 at SZAs greater 
than 78° approaching the terminator. The magnitude of the 
IBSL contamination in the radiance measurements is as high 
as 25-50% at SZA = 88°. In order to accurately retrieve 
stratospheric ozone in these situations, which frequently oc- 
cur in polar regions, it is necessary to correct the IBSL error 
in the SBUV/2 radiance measurements. This paper reports 
our analysis and characterization of the SBUV/2 IBSL error 
in orbit using N1 9 as an example, which has the largest IBSL 
error among SBUV/2 instruments for an easy illustration. 


2 SBUV/2 instrument 

Figure 3 is a schematic diagram of the SBUV/2 optical de- 
sign. SBUV/2 instruments are equipped with double Ebert- 
Fastie monochromators, which can block OBSL at an extinc- 
tion rate of 10 -6 . The monochromator covers a wavelength 


range from 160 to 405 nm, with a nominal spectral band- 
pass of 1.13 nm FWHM. SBUV/2 instruments make discrete 
measurements for profile ozone monitoring at nominal wave- 
lengths of 252.0, 273.7, 283.2, 287.7, 292.4, 297.6, 302.0, 
305.9, 312.7, 317.6, 331.3 and 339.9 nm (except for Nimbus- 
7 SBUV, which had channel 1 at 255.7 nm). It takes 24 s to 
scan through these 12 wavelengths, and completes the scan 
cycle in 32 s. The cloud cover radiometer (CCR), installed 
side-by-side with the monochromator, makes continuous and 
simultaneous measurements at 378.6 nm with a 3nm band 
pass to monitor earth surface reflectivity change during the 
spectral scan. A depolarizer in front of the monochromator 
is used to eliminate the polarization effects of atmospheric 
Rayleigh scattering, concerning strong polarity of a grating 
spectrometer that acts as an analyzer. A diffuser plate can 
be deployed to reflect the solar light into the sensor’s aper- 
ture after terminator crossing to measure the solar irradiance. 
This measurement is performed every day during normal op- 
erations. Since the radiance measurement and the irradiance 
measurement share the same monochromator and CCR, the 
sensor’s sensitivity is canceled in the albedo ratio used for 
ozone retrieval. Therefore, the bidirectional reflectance dis- 
tribution function (BRDF) calibration of the solar diffuser, 
the only different element in the optical paths between the ra- 
diance and irradiance measurements, is more important than 
the absolute photometric calibration for the ozone retrieval. 
Rigorous SBUV/2 albedo calibrations are performed in labo- 
ratory before launch, and are traceable to NIST (National In- 
stitute of Standards and Technology) BRDF standards. The 
SBUV/2 diffuser reflectivity changes, from prelaunch tests 
to on-orbit conditions, and over years of operation, are mon- 
itored using emission lines from an on-board mercury lamp. 
Typically, the diffuser reflectivity at 250 nm degraded at a 
relative annual rate of no more than 3 % per year, with sig- 
nificantly smaller changes observed in the long wavelengths. 
The diffuser calibration measurements in orbit were sched- 
uled on a weekly basis to ensure the accuracy of the albedo 
calibration. Details of SBUV/2 calibrations can be found 
elsewhere in literature (Frederick et al., 1986; Heath et al., 
1993; Hilsenrath et al., 1995; Janz et al., 1995; DeLand et 
al., 2012). 

Figure 4 shows a picture of an SBUV/2 instrument in the 
laboratory, with flight baffle plates installed surrounding the 
aperture and diffuser deployment mechanism. These plates 
are designed to shield the aperture from both direct solar 
incident radiation and spacecraft surface reflections of both 
solar light and earth shine. Most parts facing the aperture 
are anodized in black, which typically has a low reflectiv- 
ity of several percent in the UV wavelength range of interest 
(Lowery, 1977). The noticeably brighter horizontal line in 
the picture is the cylindrical axle of the diffuser deployment 
mechanism, which reflects roof light into the camera. When 
the axle is illuminated with a flood light at different angles, 
the camera can capture the light on a different part of the 
axle due to the cylindrical geometry. The axle is a couple of 
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Fig. 4. SBUV/2 instrument in the laboratory. Pieces of anodized 
plates are installed surrounding the SBUV/2 aperture as light baf- 
fles. The diffuser deployment mechanism has a cylindrical axle that 
reflects roof light into the viewer’s eye. Courtesy of Ball Aerospace 
and Technologies Corp. 


centimeters above the front baffle in the picture. In orbit, the 
aperture is oriented downward towards Earth, and the axle 
is below the front baffle. When the satellite moves from the 
daylight to the nightside near terminator, the diffuser deploy- 
ment mechanism faces the sun. At SZA = 75°, solar light can 
graze over the edges of the baffle plates, and illuminate a 
portion of the axle toward the depolarizer. The entrant slit of 
the monochromator, 1.3 mm x 30 mm in size, is behind the 
depolarizer. Light scattering and reflection from optical sur- 
faces of the depolarizer can redirect some of the light into 
the entrance slit within the field of view. This geometry is 
consistent with the finding of IBSL error onset and rapid ris- 
ing around SZA = 77° near the sunset, shown in Fig. 2. It 
also explains the clean radiance signal at sunrise when the 
diffuser deployment mechanism turns its back against the 
sun. SBUV/2 instruments use a depolarizer, which consists 
of 4 quartz wedges. In the original instrument design, the 
surfaces of the 4 quartz wedges are optically contacted with 
each other. During thermal vacuum testing of N16 in 1987, 
two of the quartz crystal surfaces separated causing a sharp 
change in the instrument’s sensitivity, which was later fixed. 
To avoid further similar problems, the last 4 SBUV/2 in- 
struments manufactured, which flew on the N14, N17, N18 
and N 1 9 spacecraft, were modified to have 1 mm air gaps 
between the surfaces of the depolarization wedges. Appar- 
ently, the air-gapped quartz wedges have surface scattering 
and reflection 1 order of magnitude larger than the optically 
contacted surfaces, which results in significant IBSL errors. 
Note that the CCR does not have the IBSL problem because 
its aperture is not covered by the depolarizer. The IBSL prob- 
lem was first uncovered in N 1 7 orbital data analysis as early 
as 2004. However, the last two SBUV/2 instruments to be 


Fig. 5. Initial profile albedo residue (IPAR) comparison between 
NOAA-19 and NOAA-17 in the Northern Hemisphere, where N19 
is contaminated with IBSL and N17 is clean. The comparison is 
made at the same latitude band of 0.5° on a weekly average. The dif- 
ference is plotted against N19’s spacecraft centered solar elevation 
angle (SCSEA). Under normal satellite attitude conditions, the SC- 
SEA is equal to SZA — 90°. Therefore, these figures cover the day- 
side from SZA = 70° to SZA = 90°. Standard errors in the average 
are plotted as error bars. IPAR differences, weighted with 1/error 2 , 
are fitted separately with two straight lines at SCSEA > — 10° and 
SCESA< —15°, shown with solid, thick black lines. The fitted 
IPAR is subtracted with the background at —15°, shown in red. 

launched on board N 1 8 and N 1 9 could not be modified be- 
cause of both flight schedules and budgetary constraints. In 
addition, there was no opportunity to perform laboratory tests 
to fully characterize the IBSL error before flight. Therefore, 
we have to use on-orbit SBUV/2 measurements to character- 
ize the IBSL error. 

3 Characterization of IBSL 
3.1 Dayside 

We have identified the IBSL error in Fig. 2 using direct com- 
parison of N16 and N18 albedo measurements at 273 nm. 
These instruments were in similar orbits in 2005, which re- 
duces possible differences in measured radiances due to di- 
urnal ozone variations. Otherwise, the measured radiances 
as a function of SZA by two instruments at different or- 
bits can be very different in their field of views, scene re- 
flectivity, terrain height and other atmospheric and geophys- 
ical properties. In the V8 SBUV ozone retrieval, an a pri- 
ori ozone profile is constructed as an initial estimate of the 
ozone profile, which has the total column ozone equal to 
the total column ozone determined with measurements at 
318 nm for ozone absorption and at 331 nm for scene re- 
flectivity (Bhartia et al., 2013). The initial ozone profile is 
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used to calculate albedo values at other short wavelengths 
in the same spectral scan, using the atmospheric backscat- 
tering forward model. The difference between the measure- 
ment and the computed albedo value, called the initial profile 
albedo residue (IPAR) in this paper, contains the difference 
between the actual ozone profile and the initial estimate, the 
temperature dependent ozone absorption cross-section error, 
the calibration offset, as well as any IBSL error. Weekly 0.5° 
latitude band zonal means of IPAR are computed and com- 
pared between two SBUV/2 instruments at the same latitude 
to minimize both geophysical and atmospheric differences. 
Figure 5 shows the weekly average differences in IPAR be- 
tween N19 and N17 near the North Pole as a function of 
the SCSEA of N19, where N1 9 approaching sunset was con- 
taminated with IBSL error and N17 at sunrise was not con- 
taminated by IBSL error. Since this error is observed near 
spacecraft sunset, N17 measurements have significant IBSL 
errors near the South Pole, while N19 measurements in that 
region are not contaminated. Therefore, a similar comparison 
to Fig. 5 near the South Pole would reveal the N17 IBSL er- 
ror. The spacecraft centered coordinates, SCSEA and space- 
craft centered solar azimuth angle (SCSAA), are preferred 
for the IBSL analysis over the reference frame on the Earth’s 
surface, since the IBSL error is attributed to the direct so- 
lar incidence on the spacecraft rather than earth shine. We 
will demonstrate the strong correlation between IBSL and 
SCSAA in the next sections, whereas there is no direct cor- 
relation between solar azimuth angle (SAA measured at the 
surface clockwise from due north) and IBSL. The IBSL er- 
ror rises rapidly in the region between SCSEA = —15° and 
SCSEA = —10° (corresponding to SZA = 75° and 80°), and 
flattens between SCSEA = —10° and SCSEA = —2° (cor- 
responding to SZA = 80° and 88°). The error bars are the 
standard errors for the latitude band zonal means. Note also 
the offset of 3 x 10 -6 at 273 nm at SCSEA = —15°, which is 
about 3 % of measured albedo value. This albedo offset can 
be a combination of calibration offsets in both N19 and N17, 
which have a 1 % uncertainty (DeLand et al., 2012), as well 
as possible diurnal changes in ozone profiles due to the dif- 
ferent local times of the measurements (Haefele et al., 2008). 
The offset fluctuated from week to week at different wave- 
lengths from —2 x 10 -6 to 7 x 10 -6 . This is larger than the 
likely calibration offset or instrument noise effect, and thus 
is likely due to atmospheric and geophysical differences. 

The IPAR differences, weighted with the reciprocal of the 
squared standard error values, are fitted with two straight 
lines separately over two SCSEA ranges, [20°,— 15°] and 
[—10°, — 8°], as shown by the thick, solid black lines. The 
fitting gives IPAR differences at SCSEA = —10° and SC- 
SEA = -15°. We estimate the IBSL error at SCSEA = -10° 
by subtracting the IPAR difference at SCSEA = —15° from 
the difference at SCSEA = — 10°, which is plotted in red. We 
also tried to use extrapolation of a straight line fitted between 
SCSEA = -20° and SCSEA = -15° for the offset subtrac- 
tion, but found that the weekly average results were too noisy 



Fig. 6. N19 IBSL error values at 273 nm at SCSEA = —10° and 
SCSEA = 6°. At SCSEA — —10°, the dayside weekly averages are 
derived from comparison of IPAR with N17. The noise (standard 
errors) in the weekly 0.5° latitude zonal means of the IPAR differ- 
ences propagated through the linear regression resulted in uncer- 
tainty in the estimated IBSL errors. Then, a piecewise linear fit is 
used to smooth and interpolate the weekly values. At SCSEA = 6°, 
N19 IBSL errors are derived directly from daily earth view mea- 
surements on the nightside (see Sect. 3.2), which have standard er- 
rors of around 0.3 %. 


to be acceptable. We do not use this IPAR difference to de- 
termine the IBSL error at SCSEA > — 10° beyond the initial 
step-up point because of the large uncertainty in the offset 
subtraction. We limit the use of the dayside data to channels 
2-4 (273.5, 283.1, 287.6 nm). Channel 1 is not included in 
this analysis because measurement noise is much more sig- 
nificant due to the low signal level at 252 nm. Determina- 
tion of IBSL error at the longer wavelengths (channels 5 and 
up) is more difficult because the measurements are impacted 
more by surface reflectivity variations and other differences 
in atmospheric properties in the lower atmosphere. The IPAR 
comparison can also be limited seasonally because of varia- 
tions in orbital geometry. The observed N1 9 IBSL error val- 
ues at SCSEA = —10° are plotted as a function of time in 
Fig. 6. The standard errors in the weekly 0.5° latitude band 
zonal means are propagated through the linear regressions, 
and produce the uncertainties for the IPAR differences. The 
estimated uncertainty in the dayside weekly IBSL errors var- 
ied from 5 to 10 %, as plotted in error bars. 

Figure 7 shows 2 yr of data points in the IPAR zonal mean 
comparison at 273 nm between N1 9 and N17. The thin solid 
curve is derived from a piece-wise linear fit (ASCSEA = 1°) 
weighted with the reciprocal of the square of the standard 
error. Note that this curve is not necessarily centered with re- 
spect to the spread of data points because of the weighting. 
A straight line is then fit to the IPAR difference between SC- 
SEA = —20° and SCSEA = —15°, and extrapolated to cor- 
rect the offset. The corrected fit result is plotted as a thick 
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Fig. 7. Two years of data points from IPAR comparisons between 
N19 and N17 in the Northern Hemisphere illustrate the IBSL er- 
ror at 273 nm as a function of SCSEA, a rapid rising edge and 
flat top before and after SCSEA = — 10°. Piece-wised linear fitting, 
weighted with statistical error of each point, is used to get a smooth 
average (thin solid curve). Then, it is corrected for background off- 
set which is an extrapolation of a straight line fitted between SC- 
SEA = —20° and SCSEA = — 15° (dashed line). The final result for 
IBSL is plotted as a thick, solid curve. 


Fig. 8. Example of N19 IBSL nightside daily measurements on the 
nightside on day 171 of 2010. Larger IBSL errors were measured 
in the Northern Hemisphere in blue, and lower signals are for the 
Southern Hemisphere, in red. Thin solid curves are for IBSL af- 
ter the empirical correction of earth shine near terminator. Thick 
straight lines are respectively fitted between SCSEA = 6° and 10° 
for daily IBSL at SCSEA = 6°, and between SCSEA =16° and 20° 
for daily IBSL at SCSEA =18°. 


solid curve. Two results come from this exercise. First, the 
IBSL error between SCSEA = -10° and SCSEA = -2° can 
be approximated well with a straight line. We will extend 
this linear approximation to SCSEA = 6° for further analy- 
sis. Second, the smooth curve between SCSEA = —15° and 
SCSEA = — 10° provides the shape of the rising edge of the 
IBSL error. 

3.2 Nightside 

As discussed above, the use of dayside data is limited to 
IBSL errors at SCSEA = — 10° (SZA = 80°) at 3 short wave- 
lengths from 273 to 289 nm, and determination of the IBSL 
shape (rapid rising edge and relatively flat top) as a function 
of SCSEA. To model the behavior of the dayside IBSL er- 
ror for all conditions, we have to look further into the night- 
side, where both the SCSAA and SCSEA dependence of the 
IBSL error can be accurately measured and characterized. 
Figure 8 shows typical nightside IBSL data measured by 
N19 in just one day. The nightside measurements are cali- 
brated in the same way as the albedo measurements in day- 
light (DeLand et al., 2012). At short wavelengths from 252 
to 298 nm, the earth shine signal drops rapidly to negligi- 
ble level at SCSEA = 6° (SZA = 96°) as measured near the 
South Pole where N19 does not have an IBSL error. At wave- 
lengths longer than 301 nm, earth shine decreases slowly to 
a negligible level at SCSEA = 9° (SZA = 99°). The mag- 
nitude of the nightside IBSL error from 252 to 340 nm in- 
creases by a factor of 2 in terms of albedo, compared to 


the larger spectral dependence of the dayside albedo. This 
weaker wavelength dependence confirms the assumption that 
the IBSL error is caused by some metal surface reflection of 
direct solar irradiance rather than earth shine. IBSL contam- 
ination in the albedo measurements at wavelengths longer 
than 318nm is negligibly small (< 0.1 %) because of the 2 
orders of magnitude increase in Earth albedo at these wave- 
lengths (Fig. 1). Similarly, there is a negligible IBSL error 
in the solar irradiance measurements, because the minimum 
solar irradiance signal at all wavelengths is at least 4 orders 
of magnitude larger than the IBSL error. 

The IBSL error is not necessarily a linear function of SC- 
SEA. As shown in Fig. 8, IBSL can be fitted well (within 3 % 
of the overall error) between SCSEA = 6° and SCSEA = 22° 
with a straight line. However, it can deviate significantly 
from extrapolations of the fitted line at SCSEA = 27°, partic- 
ularly when SCSAA is high (70° in the example). To reduce 
uncertainty, we choose IBSL measurements at SCSEA = 6°, 
as close as possible to the dayside, for interpolation of the 
IBSL error at short wavelengths. However, the full use of 
IBSL measurements at SCSEA = 6° will be limited due to 
earth shine at the long wavelengths. We will first study the 
IBSL error at SCSEA = 18°, where the measurements have 
sufficient accuracy, to understand its time dependence, go- 
niometric dependence and wavelength dependence. 
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in a reference year, defined as April 2010 to March 2011: 

Ri(t) = IBSL(C <p)/IBSL(Y r , ip) = F(t)/F(t r ), (2) 

where t r is the time in the reference year, and Eq. (1) is used 
to eliminate the goniometric factor. For a given set of t and 
<p, two consecutive days in the reference year are selected 
for having SCSAA adjacent to <p, and their IBSL errors are 
interpolated as a function of (p to obtain IBSL(/ r , cp), where t r 
is an interpolated time between these two consecutive days. 
Depending on the season, as many as three periods in the 
reference year may occur with adjacent cp values. We select 
the period with the largest time interval. Then, normalized 
drift rates are calculated, 

ARi(t)/At = (Ri(t)-l)/(t-t v ). (3) 


Fig. 9. N19 IBSL at 273 nm at SCSEA=18° is a product of 
time dependence and SCSAA dependence, (a) Daily IBSL values 
(blue dots for raw data) are highly correlated with SCSAA (dashed 
curve), and can be represented by F(t) x G(SCSAA) (solid curve), 
(b) Average drift rates are derived from comparison of daily IBSL 
values at the same SCSAA. Missing data are due to SCSAA values 
outside the range in the reference period, (c) Initial estimate of IBSL 
time dependence is derived from the comparison of daily IBSL val- 
ues at the same SCSAA, and from the average drift rate for the ref- 
erence period. It is smoothed and interpolated with piecewise linear 
fit. (d) Daily IBSL values before removing the time dependence 
are plotted against SCSAA at SCSEA =18°. (e) IBSL/F(f) is fitted 
with a polynomial of SCSAA, G(SCSAA). (f) Relative changes are 
reiterated with daily IBSL values divided by the fitted G(SCSAA), 
then normalized and smoothed for the time dependence function, 
F(t). 

3.2.1 Separation of time dependence and goniometric 
dependence 

A linear fit of daily IBSL measurements between SC- 
SEA = 14° and SCSEA = 22° gives values for both the 
IBSL error and the rate of change dIBSL/dSCSEA at SC- 
SEA = 18°. Figure 9a shows the derived daily IBSL error 
at SCSEA = 18° at 273 nm as a function of time (dots) for 
N19. The seasonal variations have a strong correlation with 
variations in SCSAA (dashed curve). Since the reflectivity of 
surfaces near the instrument entrance slit may change with 
time, and the illumination geometry of the SBUV/2 instru- 
ment will vary with season, the overall IBSL error is assumed 
to be a product of a time dependence function and a gonio- 
metric function: 

IBSL(C (p) = F(t) x G(<p), (1) 

where t is time, (p is SCSAA, and the wavelength dependence 
in F and G is implied. Note that cp is also a function of time, 
and (p, F and G are all considered at SCSEA = 18°. In order 
to derive the time-dependent factor Fit), we first compare 
each daily IBSL measurement with one at the same SCSAA 


The results are plotted in Fig. 9b. No drift rates in the ref- 
erence year are included to avoid large noise associated with 
the small time interval. There are missing days in the fig- 
ure because some SCSAA values are outside the range of the 
reference period. The average drift rate over 4 months before 
April 2010 and 4 months after March 2011 is calculated to 
determine the drift rate during the reference period, shown 
by the solid line. This allows us to make an initial estimate 
of the IBSL time dependence during the reference period, us- 
ing a linear approximation. An initial estimate of IBSL time 
dependence can then be obtained with 

fit) = [1 + it- fb) X (ARi(t)/At)] x Riit) « F(t)/F(t h ), (4) 

where (ARf it)/ At) is the average drift rate during the refer- 
ence year, % is chosen to be at the beginning of the reference 
year for convenience, and the sign could be replaced 
with the “=” sign if the drift in the reference year is strictly 
linear. Thus far, the time dependence function for IBSL is 
normalized at t\>. We renormalize it to the first day of N19’s 
earth radiance measurements, 

F n it) = fit)lfiti)KFit)IFiti), (5) 

where t\ is the time of the first day of N19’s radiance mea- 
surements. The resulting F n (t) values are shown by the dots 
in Fig. 9c. Note that F n (t) ~ F(t ) will be valid when we 
choose F(t) in Eq. (1) to be normalized on day one. These 
F n ft) values are smoothed and interpolated with a piece-wise 
linear fit to get an initial estimate of the time dependence, 
shown as the solid curve. The IBSL daily values shown in 
Fig. 9a can also be plotted as a function of SCSAA, as shown 
in Fig. 9d. Dividing these daily IBSL values by the Fit) 
function derived in Fig. 9c gives us the IBSL SCSAA de- 
pendence, which is now a well behaved function, as shown 
by the dots in Fig. 9e. A 4th order polynomial function fits 
these data well, as shown by the smooth curve, which will 
be used as Gicp) in Eq. (1). The fitted Gi<p) function is then 
used to remove the SCSAA dependence from the daily IBSL 
measurements. This creates a new estimate of the IBSL time 


www.atmos-meas-tech.net/7/267/2014/ 


Atmos. Meas. Tech., 7, 267-278, 2014 


274 


L. K. Huang et al.: Characterization of in band stray light in SBUV/2 instruments 


At SCSAA 18.0 degrees (night side) 



2009 2010 2011 2012 

Date (year) 


Fig. 10. (a) dIBSL / dSCSEA/F(t) can be fitted with a polynomial of 
SCSAA, S(SCSAA). (b) Daily values of the slope dIBSL /dSCSEA 
as a function of time (dots) can be approximated with F(t ) x 
S(SCSAAO) (solid curve). 

dependence using all measurements, as shown by the dots in 
Fig. 9f. These values are smoothed to determine the time de- 
pendence function F(t). This process is iterated once. The fi- 
nal IBSL time dependence is shown with the smooth curve in 
Fig. 9f. Combining the IBSL time dependence and the IBSL 
SCSAA dependence from Eq (1), we can calculate the IBSL 
error at SCSEA =18°, shown by the smooth curve in Fig. 9a, 
where SCSAA is measured at SCSEA =18° and is approxi- 
mately constant in a given day. 

Daily values of dIBSL /dSCSEA at 273 nm are plotted 
as a function of time in dots in Fig. 10b. We find that 
dIBSL /dSCSEA is also correlated with SCSAA. The nor- 
malized slope is shown as a function of SCSAA in Fig. 10a, 
where we have removed the time dependence F(t) shown 
in Fig. 9f. The normalized slope can be then fitted with a 
polynomial function, S(cp), shown by the smooth curve in 
Fig. 10a. Thus, we can reproduce the time dependence of the 
slope with the fitted functions, 

dIBSL/dSCSEA(0 = S(<p(t)) x F(t), (6) 

which is shown by the smooth curve in Fig. 10b. 

3.2.2 Goniometric functions at SCSEA = 18° 

All 12 wavelength channels of daily IBSL values at SC- 
SEA = 18° are processed to derive their goniometric func- 
tions in the same way as described in the previous sec- 
tion. Figure 11a shows all 12 fitted G(cp,X) functions at SC- 
SEA = 18° for N19. Since the time dependence functions, 
F(t, X), are normalized on the first day of N19’s earth 
view radiance measurements, G(cp, X) gives IBSL at SC- 
SEA = 18° on that date when cp = 6 1.1°. The fitted polyno- 
mials in Fig. 1 la are very similar to each other in shape. They 



Spacecraft-Centered Solar Azimuth Angle At SCSEA 18 degrees (degrees) 


N19 IBSL SCSAA Dependence Scale Factor 



wavelength (nm) 

Fig. 11. Dashed curves in the top panel are for 12 channels of fit- 
ted G(SCSAA) functions, which can be approximated within 0.7 % 
of IBSL by a product of a wavelength independent factor for the 
SCSAA dependence and the scale factors in the bottom panel, 
go(SCSAA) xC(X). 



Spacecraft-Centered Solar Azimuth Angle At SCSEA 1 8 degrees (degrees) 

Fig. 12. Fitted S(SCSAA) functions for the IBSL slope, 
dIBSL /dSCSEA, at different wavelengths can be approximated 
with an average (solid line). 

can be approximated by 

g(cp,X) = g 0 ((p) x C(X), (7) 

where go((p) is an average of G(cp, X) over all wavelengths, 
shown as the thick solid curve, and SCSAA is taken at SC- 
SEA = 18°. The scale factors C( X), shown in Fig. lib, then 
minimize the standard deviation between the scaled observa- 
tions and the average. The differences between G(cp, X) and 
g(cp, X) are less than 0.7 % of the IBSL error, except for chan- 
nel 1 with a maximum of 1 .2 % due to low signal level at 
252 nm. Therefore, g(cp, X) is a valid approximation of the 
IBSL error’s goniometric dependence at SCSEA = 18°. 
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Figure 12 shows the fitted S(<p) functions for the slope 
dIBSL/dSCSEA normalized by the time dependence F(t) 
at all 12 wavelengths. We can see that S(cp) is essentially 
wavelength-independent. An average of S(cp) is computed 
(thick solid curve). A maximum difference of 10 -7 at SC- 
SAA = 73° is found between the average and channel 12, 
which is less than 0.5 % of the IBSL error. Therefore, an av- 
erage of the fitted S (cp) functions can be used for the gonio- 
metric dependence of the slope at all wavelengths. 

From the above exercises we conclude that the goniomet- 
ric functions for the IBSL and dIBSL/dSCSEA can be ap- 
proximated as the combination of a wavelength-independent 
angular function and a wavelength-dependent scale factor. 
This is because the light baffle shape determines how the 
baffle blocks the solar light in various solar incidence an- 
gles, which is purely geometrical. The reflection of illumi- 
nated surfaces can also have both the wavelength dependence 
and the angular dependence. The separation of the wave- 
length dependence and the angular dependence, as described 
above, provides a convenient and suitable approximation for 
the practical purpose of IBSL correction. 

3.2.3 IBSL at SCSEA = 6° 

The radiance measurements before the sunrise in the South- 
ern Hemisphere, by N19 near the South Pole as the exam- 
ple in Fig. 8, show how ozone absorption blocks the solar 
light in the UV spectrum. At SZA = 99°, the albedo values at 
340 nm in Southern Hemisphere is about 1 % of IBSL found 
in the Northern Hemisphere, where the earth shine is from 
light scattering of atmosphere above 90 km. It dropped below 
the sensor’s noise level at SZA= 100° in Southern Hemi- 
sphere, which corresponds to a height of 100 km. (N19 picks 
up some negligible IBSL in Southern Hemisphere starting 
at SCSEA = 18°.) Ozone absorption cross sections at wave- 
lengths shorter than 300 nm are 2-3 orders of magnitude 
larger than that at 340 nm. Therefore, solar light at the short 
wavelengths are completely blocked by ozone, rather than 
by an Earth eclipse, at much lower altitude before reach- 
ing SZA = 96°. To estimate IBSL around SCSEA = 6°, an 
empirical correction of the earth shine after the sunset is at- 
tempted, which scales the earth shine measurements before 
the sunrise with CCR values at SZA = 96°: 

7sunset(V SZA) — Cunris e(^> SZA) 

X [/sunset(A = 378 nm, SZA = 96°)/ 
f sunrise = 378 nm, SZA = 96°)], (8) 

where 7 SU nS et and / SU nrise are in the albedo values. This is con- 
sidered as a zero-order approximation. The resulting IBSL 
after the earth shine correction are plotted in thick solid 
curves in Fig. 8. This earth shine correction has no impact 
at the short wavelengths at SCSEA = 6°. At the long wave- 
lengths, residues of earth shine correction contribute to noise 
in daily IBSL values at SCSEA = 6°, about 1-2% in terms 


Time Dependence Function F(t) for N1 9 IBSL at SCSEA 6° 



Fig. 13. Time dependence functions F(t,X ) for N19 IBSL at SC- 
SEA = 6°, normalized to the start of the mission. The magnitude 
of the overall drift monotonically increases towards shorter wave- 
lengths. At the end of the data record, the magnitude of F(t,X) 
ranges monotonically from 0.68 at 252 nm to 0.90 at 340 nm. 


of standard deviation (with a peak-to-peak noise of 10% of 
IBSL), which is negligible in comparison with albedo values 
in daylight. The daily IBSL values at SCSEA = 6° are de- 
rived from a smooth function fit over SCSEA = 6° and SC- 
SEA = 12°. The resulting daily IBSL values at 273 nm are 
plotted in Fig. 6. 

Estimated daily IBSL values at SCSEA = 6° are then pro- 
cessed for the separation of time dependence and gonio- 
metric dependence terms in the same way as described in 
Sect. 3.2.1. Figure 13 shows the time dependence functions, 
F(t,X), at all 12 wavelengths. The results at SCSEA = 6° 
are consistent with those at SCSEA = 18°. For example, at 
273 nm the calculated relative IBSL error at SCSEA = 6° de- 
creased to 0.751 at the end of August 2011, compared to a 
calculated value of 0.775 at SCSEA = 18° shown in Fig. 9f. 
Note that long wavelength channels have a lower drift rate 
than the short wavelength channels. This is a typical behavior 
in metal surface reflectivity in the UV spectrum in the pro- 
cesses of photo carbonization of prelaunch contamination on 
the spacecraft (Hilsenrath et al., 1994). 

The calculated SCSAA dependence of the IBSL error as 
a function of wavelength at SCSEA = 6°, G(cp, A.), is plotted 
in Fig. 14a. These curves share a similarity in shape con- 
sistent with Fig. 11a. However, only the short wavelength 
channels 2-6 are close to each other in shape to a high de- 
gree of accuracy (0.7 %). Large differences are found with 
the other channels due to the limitation at SCSEA = 6°, earth 
shine contamination at long wavelengths and low solar flux at 
252 nm. However, the percentage of IBSL errors at the long 
wavelengths is negligibly small, and the 252 nm channel is 
not used in the V8. 6 SBUV/2 profile retrieval (Bhartia et al., 
2013). Therefore, we take the average only over channels 2-6 
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Fig. 14. (a) Dashed curves in the top panel are for 12 channels of fit- 
ted G(SCSAA) functions, which monotonically increases as wave- 
length increases. Only channels 2-6 are averaged for the wavelength 
independent factor go(SCSAA) because of the earth shine contam- 
ination at the long wavelengths and low signal at 252 nm. (b) The 
scale factors in the bottom panel, C(X), result from the least square 
between G(SCSAA) and g 0 (SCSAA) xC(A). 


for the wavelength independent shape factor, go 0/0, shown 
with a thick, solid curve. The scale factors, C( X) shown 
in Fig. 14b, result in least squares of deviation between 
go(cp) x C( X) and G(cp, X), which are less than 0.25% in 
terms of standard deviations for channels 2-6. Selection of 
these 5 short wavelength channels for the average guarantee 
the accuracy for IBSL error corrections in the ozone profile 
retrieval. We believe this average is a better representation 
of the goniometric dependence for the other long wavelength 
channels and that at 252 nm. The wavelength dependence of 
C (A.) at SCSEA = 6° is characteristically comparable to that 
at SCSEA =18° shown in Fig. 11. 

The estimation of the IBSL error in daylight (SZA <88°) 
has to be based on interpolation of the derived values at 
SCSEA = -10° and SCSEA = 6°. The weekly IPAR results 
at SCSEA =—10° are smoothed and interpolated to daily 
values with the piecewise linear fitting function, as shown 
in Fig. 6 for N19 at 273 nm. The slope, dIBSL/dSCSEA, 
is then computed with straight lines connecting the two 
points. Next, the time-independent normalized slope at SC- 
SEA = 6°, S(<p,k) =(dIBSL/dSCSEA)/F(f), is calculated, 
as shown in Fig. 15 for N19 as an example. As discussed 
in Sect. 3.1, we can only obtain valid dayside IBSL at the 
6 short wavelengths. Among them, channel 1 has low sig- 
nal, and channels 5 and 6 have some scene noise. Therefore, 
relatively more reliable results of S(cp) can be obtained only 
with channels 2-4. The difference in S(cp) among channels 
2-4 equates to only 2 % of the IBSL error at the midpoint 
between SCSEA =—10° and SCSEA = 6°. Therefore, the 
wavelength dependence of S(cp) between SCSEA =—10° 


and SCSEA = 6° is negligible at least among these short 
wavelength channels. This is consistent with the finding at 
SCSEA = 18° that the normalized slopes are approximately 
wavelength independent. Therefore, the average of S (cp) over 
channels 2—4 is calculated and used to represent the normal- 
ized slope at all wavelengths. 

3.3 IBSL corrections 

Based on the discussions and results in the previous sec- 
tions, we can model IBSL between SCSEA = —10° and SC- 
SEA = 6° with the following function: 

IBSL(SCSEA, t, A.) = 

[go(cp(t)) x C(X) + S(cp(t)) x (SCSEA - 6)] x F(t, X), (9) 

where the goniometric functions, go , C and S, and the time 
dependence function, F , are derived and defined at SC- 
SEA =6°, and SCSAA = (p{t) is measured at SCSEA = 6° 
and is approximately constant in a given day. An average of 
the rising edge of the IBSL error between SCSEA = —15° 
and SCSEA =—10° over multiple years, which is derived 
as described in Sect. 3.1, is scaled to be connected with the 
straight line of IBSL(SCSEA,f,A.) at SCSEA = -10°. 

We have examined data from all SBUV/2 instruments for 
IBSL errors. N19 leads (a peak IBSL of 2.4 x 10 -5 at 273 nm 
in terms of albedo), followed by N17 (1.9 x 10 -5 ), N14 
(1.6 x 10 -5 ) andN18 (1.2 x 10 -5 ). The instruments that use 
a depolarizer with optically contacted prism elements (N9, 
N1 1, N16) have IBSL errors about 20 times less than the in- 
struments that use a depolarizer with air-gapped prisms (N14, 
N17, N18, N19). Therefore, the IBSL contamination in N9, 
Nil and N16 albedo measurements is less than 1 % even at 
SZA = 88°, which is not significant for practical purposes. 
We have used the model described here to develop sepa- 
rate empirical IBSL correction functions for N14, N17, N18 
and N19. These corrections are implemented in the ozone 
processing code for the current NOAA SBUV/2 operational 
data processing, as well as for the SBUV/2 V8.6 reprocess- 
ing product (DeLand et al., 2012, Bhartia et al., 2013). The 
estimated uncertainty of the IBSL correction is less than 
5 % of the error between SZA = 80° and SZA = 88°, which 
represents about 1 % in the earth radiance measurements at 
273 nm. 

Figure 16 shows an example ofN19 albedo measurements 
at 273 nm before and after IBSL correction. The data were 
taken on the spring equinox 2010 near the North Pole. While 
the IBSL error was approximately constant in absolute terms 
(— 2 x 10 -5 between SZA = 80° and 88°), the percentage 
contamination in the albedo measurements increases from 
-25% at SZA = 80° to -40% at SZA = 88°. Figure 17 
compares retrieved ozone profiles before and after the IBSL 
correction. The impact of uncorrected IBSL error to the re- 
trieved ozone profile above 1 hPa (Layer 1 6) is more than 
30 %. While the retrieved total column ozone amount, which 
depends on the long wavelength channels, is not affected 
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Fig. 15. N19’s dIBSL/dSCSEA at SCSEA = 6° is calculated with 
IBSL at SCSEA = —10° in the IPAR analysis on the dayside and 
measurements at SCSEA = 6° on the nightside. 
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Fig. 16. Sample of N19 albedo measurements at 273 nm near the 
North Pole before and after the IBSL correction. 

by IBSL, the V8.6 ozone retrieval algorithm pushes up the 
ozone layers below lOhPa by 1-2% to compensate for the 
deficits in the upper layers. Even though the IBSL error im- 
pacts only a small portion of the SBUV/2 data at very high 
SZA, this error can result in large errors in the retrieved 
stratospheric ozone values in polar regions. The methodol- 
ogy developed in this paper reduces the IBSL-related uncer- 
tainty in the retrieved upper layer ozone to only a few per- 
cent. 


Fig. 17. The sample of IBSL effects in the retrieved ozone profiles 
near the North Pole is taken from an average of two weeks around 
the spring equinox 2010. Twenty-one SBUV V8 profile layers are 
defined by 10~ 2w atmosphere pressure levels, where n starts from 0 
for the bottom boundaries. Only top layers above 1 0 hPa are labeled. 
Other lower layers overlap each other within 2 %. The black curve 
for the total column ozone difference is near zero. 

4 Conclusions 

IBSL is found in the earth radiance measurements from the 
last 4 SBUV/2 instruments, which were flown on the NOAA- 
14, 17, 18 and 19 satellites, at a SZA larger than 77° when 
approaching terminator at sunset. The source of IBSL is from 
direct solar illumination of the spacecraft rather than from 
earth shine. Large surface reflection and scattering that are 
associated with the air-gapped depolarizer results in signifi- 
cant IBSL errors. We have characterized the IBSL error be- 
tween SCSEA =—10° and SCSEA = 6° (corresponding to 
SZA = 80° and 96°) with a linear function of SCSEA. The 
magnitude of IBSL is characterized as a function of both 
the SCSAA and the mission time. We have implemented an 
IBSL correction in the V8.6 SBUV/2 ozone retrieval with an 
uncertainty of 1 % of albedo values. The maximum uncor- 
rected IBSL error can reach 40 % of albedo values at 273 nm 
at SZA = 88°. The retrieved layer ozone above 1 hPa can 
be underestimated by more than 30%. The IBSL error cor- 
rection reduces the IBSL-related uncertainty in the retrieved 
stratosphere ozone profile to a few percent. 
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